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We show that the naked singularities arising in dust collapse from smooth initial data (which 
include those discovered by Eardley and Smarr, Christodoulou, and Newman) are removed when 
we make transition to higher dimensional spacetimes. The cosmic censorship is then restored for 
dust collapse which will always produce a black hole as the collapse end state for dimensions D > 6, 
under conditions to be motivated physically such as the smoothness of initial data from which the 
collapse develops. 
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A considerable debate has continued in recent years on 
the validity or otherwise of the cosmic censorship conjec- 
ture (CCC) introduced by Penrose and Hawking, which 
effectively states that any singularities that arise from 
gravitational collapse from a regular initial data must 
not be visible to far away observers in the spacetime and 
are always hidden within black holes. Such an assump- 
tion has been used extensively and is fundamental to the 
theory as well as applications of black hole physics. On 
the other hand, if a naked singularity results as collapse 
end state, it is no longer necessarily covered by the event 
horizon and could communicate, in principle, with out- 
side observers. Such a scenario is of physical interest be- 
cause a naked singularity may have theoretical and obser- 
vational properties quite different from a black hole end 
state, and communications from extreme strong gravity 
regions dominated by quantum gravity may be possible. 

Though there is no satisfactory proof or mathemati- 
cal formulation of CCC available despite many efforts, 
there are many examples of dynamical collapse models 
available which lead to a black hole (BH) or a naked sin- 
gularity (NS) as the collapse end state, depending on 
the nature of the initial data (see e.g. Q and refer- 
ences therein). In particular, pioneering analytic models 
by Christodoulou, and Newman, and numerical work by 
Eardley and Smarr , established the existence of shell- 
focusing naked singularities as the end state of a continual 
collapse, where the physical radius of all collapsing shells 
vanishes. In these models the matter form was taken 
to be marginally bound dust, assuming that the initial 
data functions are smooth and even profiles. Newman 
generalized these models for non-marginally bound class, 
thus covering the entire class of dust collapse solutions. 
A general treatment for dust collapse with generic initial 
data was developed in Q . 

In fact, gravitational collapse of a dust cloud has 
served, over past several decades, as a basic and fun- 
damental paradigm in black hole physics. This was first 
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discussed by Oppenheimer and Snyder j^j, who consid- 
ered a continual collapse of a homogeneous spherical dust 
cloud in general relativity. The important feature that 
emerges is that as the collapse progresses, an event hori- 
zon of gravity forms, which fully covers the spacetime sin- 
gularity developing later as the collapse end state. Since 
the singularity is always covered by the horizon, we say a 
black hole has developed as the collapse final state. This 
scenario was the basic motivation and origin for CCC, 
stating that any physically realistic collapse must have a 
similar outcome producing a black hole only. 

Our purpose here is to show explicitly, that the above 
mentioned naked singularities of dust collapse (for both 
marginally bound and non-marginally bound dust) are 
removed when we go to higher (> 6) dimensions, if we 
allow for only smooth and analytic initial profiles as as- 
sumed by |2| • In other words, the collapse of a dust- 
like matter, whose initial density and velocity profiles are 
taken to be smooth C°° functions, which is a condition 
to be motivated physically, always goes to a black hole 
end state for dimensions D > 6, and CCC is restored. 
We point out that the increase in dimensions deforms 
the apparent horizon, which is the outer boundary of 
trapped surfaces, in such a manner that the formation of 
the trapped surfaces is advanced. Hence, there exists a 
critical dimension beyond which the neighborhood of the 
center gets trapped before the central singularity and the 
final outcome of collapse is a black hole necessarily. 

Such a scenario offers a possibility to recover CCC, be- 
cause it is quite possible that we may be actually living 
in a universe with higher dimensions. The recent devel- 
opments in string theory and other field theories strongly 
indicate that gravity is possibly a higher dimensional in- 
teraction, which reduces to general relativistic descrip- 
tion at lower energies. Hence, it could be that CCC is 
restored due to extra physical effects arising from our 
transition itself to a higher-dimensional spacetime con- 
tinuum. The recent revival of interest in this problem is 
motivated partly by the Randall-Sundrum brane-world 
models |f|. 

To study the collapse of a spherical dust cloud, we 
choose a general spherically symmetric comoving metric 
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in N > 4 dimensions which has the form, 

ds 2 = -e^dt 2 + e 2 ^dr 2 + R 2 (t, r)dQ 2 N _ 2 (1) 
where, 
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is the metric on (TV — 2) sphere. The energy-momentum 
tensor of the dust has the form, 
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We take the matter field to satisfy the weak energy con- 
dition, that is, the energy density measured by any local 
observer be non-negative, and so for any timelike vector 
V 1 , we must have, 
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In the case of a finite collapsing cloud, there is a finite 
boundary < r < r&, outside which it is matched to a 
Schwarzschild exterior. The range of the coordinates for 
the metric is then < r < oo, and — oo < t < t s (r) 
where t s (r) corresponds to the singular epoch R = 0. 
Solving the iV-dimensional Einstein equations we get the 
generalized Tolman-Bondi-Lemaitre (TBL) metric @ as, 
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dr 2 + R 2 (t,r)dfl 2 N _ 2 (5) 



where f(r) is an arbitrary function of the comoving ra- 
dius r, and f(r) > —1. Also we get the required equa- 
tions of motion as, 
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Here F(r) is another arbitrary function of the comoving 
coordinate r. In spherically symmetric spacetimes F(r) 
has the interpretation of mass function which describes 
the mass distribution of the dust cloud and f(r) is the 
energy function, which specifies the velocity profiles for 
the collapsing shells for any given mass function. The 
energy condition then implies F 1 > 0. It follows from the 
above that there is a spacetime singularity at R = and 
at R' = 0. The latter are called 'shell-crossing' singular- 
ities, which occur when successive shells of matter cross 
each other. These have not been considered generally 
to be genuine spacetime singularities, and possible ex- 
tensions of spacetime have been investigated through the 
same @- On the other hand, the singularity at R = is 
where all matter shells collapse to a zero physical radius, 
and hence this has been known as a 'shell-focusing' sin- 
gularity. The nature of this singularity has been investi- 
gated extensively in four-dimensional spacetimes (see e.g. 
references in pj), and it is known for the case of spheri- 
cal dust collapse that this can be both naked or covered 



depending on the nature of the initial data from which 
the collapse develops. The collapse condition is R < 0. 

We now use the scaling independence of the comoving 
coordinate r to write (see e.g. Q), 
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This means we have scaled the coordinate r in such a 
way that at the initial epoch R = r, and at the singu- 
larity R = 0. It should be noted that we have R = 
both at the regular center r = of the cloud, and at the 
spacetime singularity, where all matter shells collapse to 
a zero physical radius. The regular center is then distin- 
guished from the singularity by a suitable behaviour of 
the mass function F(r) so that the density remains finite 
and regular there at all times till the singular epoch. The 
introduction of the parameter v as above then allows us 
to distinguish the spacetime singularity from the regular 
center, with v — 1 at the initial epoch, including the cen- 
ter r — 0, which then decreases monotonically with time 
as collapse progresses to value v = at the singularity 
R = 0. 

In order to ensure the regularity of the initial data, it 
is evident from the equations of motion that at the initial 
epoch the two free functions F(r) and f(r) must have the 
following forms, 



F{r) =r { - N -^M{r); f{r)=r 2 b(r) 



(9) 



Following , we now assume that the initial density and 
energy functions p{r) and fir) are smooth and even, en- 
suring their analytic nature. We note that the Einstein 
equations as such do not impose any such restriction, 
which has to be physically motivated, and it implies a cer- 
tain mathematical simplicity in arguments to deal with a 
dynamical collapse. It follows that both A4(r) and b(r) 
are now smooth C°° functions, which means the Taylor 
expansions of these functions around the center must be 
of the following form, 



and 



M(r) =M a + r 2 M 2 + r 4 M 4 + 



b(r) =b + r 2 b 2 + r 4 b 4 + 



(10) 



(11) 



that is, all odd terms in r vanish in these expansions, and 
the presence of only even terms would ensure smoothness. 

To predict the final state of collapse for a given initial 
mass and velocity distribution, we shall study below the 
singularity curve resulting from the collapse of successive 
matter shells, and the apparent horizon developing in the 
spacetime. A decreasing apparent horizon in (t, r) plane 
then is a sufficient condition for a BH as it shows the 
entrapment of the neighborhood of the center before the 
central singularity. 
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Consider now the endless collapse of a dust cloud to a 
final shell- focusing singularity at R = 0, where all matter 
shells collapse to a zero physical radius. With the regular 
initial conditions as above, equation© can be written as, 
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Here the negative sign implies that v < 0, that is, the 
matter cloud is collapsing. Integrating the above equa- 
tion with respect to v, we get, 
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We note that the coordinate r is to be treated as a con- 
stant in the above equation. Expanding t(v, r) around 
the center, we get, 



t(v, r) = t(v, 0) + rX(v) + r 2 ^ + r 3 ^ 



(14) 



Now from equation^UJ) and (|llf) we have X{v) = 0. The 
function X 2 (v) is then given as, 
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where, 

b = b(0);M o = M(0);b 2 = b"(0);M 2 = M"(0) (16) 

We note that the value of X 2 is determined fully by the 
initial values of the mass function M (r) and the energy 
function b(r). Thus, the time taken for the central shell 
to reach the singularity is given by 
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From the above equation it is clear that for t So to be 
defined, 



Mo + v^ N ~ 3 h > 
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that is, the continual collapse condition implies the posi- 
tivity of the above term. Hence the time taken for other 
shells close to the center to reach the singularity can now 
be given by, 



t s (r) = t So + r 



2 X 2 (0) 
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In order to determine the visibility or otherwise of the 
singularity at R = 0, we need to analyze the causal struc- 
ture of the trapped surfaces and the nature and behaviour 
of null geodesies in the vicinity of the same. If there exist 
future directed null geodesies with past end point at the 
singularity, which go out to faraway observers in space- 
time, then the singularity is naked. In the case otherwise, 



we have a black hole resulting as the end state of a con- 
tinual collapse. The boundary of the trapped region of 
the space-time is given by the apparent horizon within 
the collapsing cloud, which is given by the equation, 
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Broadly, it can be stated that if the neighborhood of the 
center gets trapped earlier than the singularity, then it 
is covered, otherwise it is naked with families of non- 
spacelike future directed trajectories escaping away from 
it. For example, it follows from the above equation that 
along the singularity curve t — t s (r) (which corresponds 
to R = 0), for any r > we have F(r) going to a constant 
positive value, whereas R — ► 0. Hence it follows that 
trapping already occurs before the singularity develops 
at any r > along the singularity curve t s (r). 

What we need to determine now is when there will be 
families of non-spacelike paths coming out of the central 
singularity at r = 0, t = t s (0), reaching outside observers, 
and when there will be none. The visibility or otherwise 
of the singularity is decided accordingly. By determin- 
ing the nature of the singularity curve, and its relation 
to the initial data, we are able to deduce whether the 
trapped surface formation in collapse takes place before 
or after the central singularity. It is this causal structure 
that determines the possible emergence or otherwise of 
non-spacelike paths from the singularity, and settles the 
final outcome in terms of either a BH or NS. From equa- 
tions (J20J), we have, 



Vah{r) = [r 2 X(r)]»-3 
Using the above equation in i|13fl we get 
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As we are considering the behaviour of the apparent hori- 
zon close to the central singularity at r = 0, R = (all 
other points r > on the singularity curve are already 
covered, as we pointed out above), therefore the upper 
limit of integration in the above equation is small, and 
hence we can expand the integrand in a power series in 
v, and keep only the leading order term, which amounts 
to, 
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Since the apparent horizon is a well-behaved surface for 
a spherical dust collapse, hence we can say that the sin- 
gularity curve for the collapse and its derivatives around 
the center are also well-defined, as the same coefficients 
are present in both (|2*3l and l(T4l . 

It is now possible to analyze the effect of the num- 
ber of dimensions on the nature and shape of the appar- 
ent horizon. As we increase the number of dimensions 
and go to dimensions higher than five, the negative term 
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FIG. 1: The apparent horizon in different spacetime dimen- 
sions. Here M{r) = 1 - 3.6r 2 H and b(r) = 1 - 3.6r 2 H , 

and apparent horizon curves are given for dimensions 4 to 7. 

in equation starts dominating, thus advancing the 
trapped surface formation in time. Fig. 1 shows the 
behaviour of the apparent horizon curve for the same 
initial data in different dimensions. As seen here, in 
usual four-dimensional spacetime the above initial pro- 
file ensures an increasing apparent horizon and it can 
be explicitly shown that the end state is a NS. But in 
six and higher spacetime dimensions we see the negative 
term starts dominating, purely as a result of increase in 
dimensions. This causes the trapped surfaces to form 
sufficiently early to cover the singularity. It follows that 
the number of spacetime dimensions causes an intriguing 
effect on the causal structure and nature of the apparent 
horizon, as seen above clearly. 

We thus see that for a smooth initial data and for di- 
mensions higher than five, the apparent horizon becomes 
a decreasing function of r near the center. This implies 
that the neighborhood of the center gets trapped before 
the central singularity and the central singularity is then 



always covered, as opposed to the four-dimensional mod- 
els by Eardley and Smarr, Christodoulou and Newman, 
in which case a positive ^(0) would make the apparent 
horizon curve increasing and hence ensures a naked sin- 
gularity. To be specific, suppose there is a future directed 
outgoing null geodesic coming out from the central sin- 
gularity at R = 0,r = 0. If (ti,n) is an event along the 
same, then t\ > t So and r\ > 0. But for any such r\, 
the trapped region already starts before t = t So , hence 
the event (ii , r\ ) is already in the trapped region and the 
geodesic cannot be outgoing. Thus, there are no out- 
going paths from the central singularity, making it cov- 
ered. It now follows in general (i.e. for both marginally 
bound and non-marginal cases) that for a dust collapse 
with smooth initial profiles, the final outcome is always a 
black hole for any spacetime dimensions D > 6, and all 
the naked singularities occurring in lower dimensions are 
removed. 

In five dimensions we have an interesting scenario aris- 
ing (see also 0). As it is clearly seen from the equa- 
tion i|23|) . we have a critical value of <%2(0), below which 
the apparent horizon is decreasing and we will get a BH 
end state. However, in the case otherwise, a naked sin- 
gularity can result. It is interesting to note also that 
the above results holds if the initial profiles, instead of 
being absolutely smooth C°° functions, are taken to be 
sufficiently smooth (i.e. at least a C 2 function). 

We note that it may still be possible to have both 
BH/NS outcomes as collapse end states when initial data 
and metric are not assumed necessarily to be C°° or at 
least C 2 analytic functions 0. Hence such smoothness 
conditions need to be investigated and probed further 
carefully to see if they could be strongly motivated from 
a physical perspective. It will also be interesting to see 
the effects of introducing pressures within the collaps- 
ing cloud, and to examine if we could still avoid naked 
singularities and preserve cosmic censorship in higher di- 
mensions by introducing physically motivated conditions 
such as above. This will be discussed elsewhere. 
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